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For an fee crystal with central force interactions and separately for a 

scalar model on a square lattice , we compute exactly the phonon transmission 

coefficient T{oj) through a disordered planar interface between two identical 

semi - infinite leads. At high frequencies T(u) exhibits a strong frequency 

dependence which is determined by the correlation length of the disorder. 



The problem of energy transfer between two solids occurs in a wide range of physical 
structures, including semiconductor quantum wells and superlattices. Early studies have 
assumed a perfect solid-solid interface and as a first approximation, an elastic continuum 
model accounting for acoustic mismatch has been employed to quantify heat flow through the 
interface |TJ]. Such a continuum theory describes generic low-frequency properties, but e.g. 
to account for experimental observations of the Kapitza thermal conductance between solids 
at high temperatures, a lattice-dynamical approach || is required, which incorporates 
important parameters such as the phonon dispersion. 

The theory of reference [§] provides a comprehensive lattice-dynamical description of 
perfect solid-solid interface. In this paper we discuss one extra mechanism which might play 
a significant role in energy transfer, namely phonon scattering due to interfacial disorder. 
Breaking the translational invariance at the interface gives rise to the so-called scattering- 
mediated phonon transmission(reflection) [Q , for which calculations based on pertubation 
theory predict a strong frequency dependence with a power-law crossover determined by the 
disorder characteristics [||. We study these effects by analysing a lattice-dynamical model 
and calculate exactly the phonon transmission coefficient as a function of frequency. The 
disorder is introduced as a correlated or uncorrelated random variation of masses m- } along 
the lattice sites j of the interface layer. 

We formulate the problem of phonon transmission across a disordered solid-solid inter- 
face by considering two identical semi- infinite leads attached to a 'scattering region'. The 
two semi-infinite perfect crystals are envisaged as waveguides for incident and scattered 
phonons. The scattering region(disordered interface) consists of a single atomic plane where 
interactions with the impurities take place. A typical geometry is shown in Fig.[|. 
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FIG. 1. A typical scattering geometry, where L(R) refer to sites on the faces of left(right) lead. 

As a model we take an fee lattice of masses m- 3 coupled to each other by central force 
springs of stiffness k. As a second example we also consider a square lattice described by a 
Born model with a single degree of freedom per site. In the latter, phonon polarization is 
neglected. In the harmonic approximation the lattice-dynamics in both cases is described 
by the linearized equations of motion of atomic displacements Uj 

mjii? = (!) 

where K?f = ft(£y£y — 4<5ij<5 a/3 ), iy = (j — i)°/|j — i|, and (j — i) is taken from the first 
coordination sphere. The indices a,/3 run through x,y,z when phonons can be polarized 
and only through z when they are treated as scalars. 

In the semi-infinite leads all atoms are of unit mass, whereas the masses at the interface 
have mean value (rrij) = 1 and variance a ||. To generate an ensemble of N positive masses, 
we first introduce the following random numbers Xj = J2k{ a k cos(kj) + bf~ sm(kj)) where 
and bk are Gaussian random numbers with zero mean value and (aka^) = (bkbk>) = &kw e ~ 
and j labels a mass along the interface. Starting from these correlated quantities with 
correlation length £, which may assume positive and negative values, we define the non- 
negative quantities rhj and rhj via the relations: rhj = Xj~X 1 j im i — a y^sy^p^- Finally 
the set of correlated mass {rrij} used in the simulation is defined by rrij = 1 — rhj + (rhj). 
In what follows, to isolate the effect of correlations, we compare the transmission coefficient 
of such an interface, with that of an uncorrelated interface obtained by randomly 'shuffling' 
the above set. 

To calculate the overall phonon transmittance we compute the unitary scattering matrix 



S(lv) at a fixed frequency using the Landauer-Biittiker formalism . In addition to various 
tests on the numerical code, such as the calculation of the density of states or of isotropic 
scattering for a single mass defect, the unitarity of S(u) is checked in all calculations. 
SS* = 1 reflects the conservation of flux. The S-matrix is extracted from the Green's 
function(Fig.[I|) 
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yG(R,L;u) G{R,R;u) , 

G{u) = {l-g{u)D.)- l g{u), (3) 

where g(u) is the Green's function for the leads calculated using an algorithm developed 
in [§. D s is the dynamical matrix D = M~ l / 2 KM- 1 / 2 , M~f = m^tfy, which includes 
all the couplings between atoms on the interface and between the scattering region and the 
leads. Eq.|3| is Dyson's equation written in a convenient form. 

In Fig.^|(a),(b) the frequency dependence of the overall transmission coefficient T(u) 
averaged over 10 disorder realizations of the interface masses for the fee vectorial Eq.|l| and 
over 50 for the square lattice is shown. The lateral width of the interface is 29 and 100 
in units of the nearest-neighbour spacing. The reason that we used the simplified scalar 
model is apparent. The computational time is quite large if we include phonon polarization 
because of the increase in the degrees of freedom. This prevents us from considering wider 
structures and eliminiting finite-size effects evident in the stucture of the plotted normalised 
to unity transmittance in Fig.0(a). In Fig.^(c) the qualitative behaviour of T(u) based on 
the results of pertubation theory || for both correlated and uncorrelated non-ideal interface 
layer is shown. These estimates are in good qualitative agreement with the results of the 
lattice-dynamical calculations. 
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FIG. 2. Phonon transmission coefficient for the (a)fcc model(£ = 5.65), (b)square scalar 
model(£ = 20), (c)estimates based on pertubation theory. In the (a) and (b) insets the corre- 
sponding normalised DOS versus uj is plotted. Frequency is measured in units of sjn/m. Inset (c) 
shows a single realization of rrij versus j for a line of correlated and uncorrelated masses (a = 0.3 
in all plots). 

Figure 2 shows that a disordered interface of just a single atomic plane gives rise to 
strong scattering-controlled phonon transmission. The overall phonon transmission coeffi- 
cient (reflection, R{uS) = 1 — T{uj)) exhibits strong frequency dependence with increasing 
incident phonon frequency. We also clearly demonstrate that the precise frequency depen- 
dence is determined by the characteristics of disorder. In particular, there exists a pro- 
nounced difference between the correlated and the uncorrelated disorder configuration. The 
exact T{u) exhibits much slower frequency dependence for correlated distribution of masses 
rrij on the plane of the solid-solid interface than for uncorrelated configuration obtained by 
shuffling the same set {rrij}. Such an effect originates from the restrictions in the phase 
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volume available for the scattered states due to the correlation-induced finite width of the 
disordered spectral distribution. For the future it would be interesting to simulate a finite 
thickness interface between dissimilar solids. 
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